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CLOSED-FORM APPROXIMATE SOLUTIONS TO THE OSCILLATORY 
PLASMA CONTINUITY EQUATIONS 
by Lester D. Nichols and J. Reece Roth 
Lewis Research Center 

SUMMARY 

Recent experimental observation of continuity- equation plasma oscillations has made 
desirable an improved closed-form approximation to the nonlinear differential equations 
which describe them. Such an approximation has been found, and is compared with the 
less satisfactory approximations previously obtained, and also with the exact solutions 
to these equations obtained by numerical integration. 


INTRODUCTION 


When the mean free paths of electrons and neutrals in a partially ionized gas are 
comparable to or larger than its dimensions, Fick's law of diffusion is inappropriate. 
One therefore cannot approximate the divergence terms of the continuity equations by the 
product of a diffusion coefficient and the Laplacian of the number density. In this case, 
the continuity equations assume the general form discussed in references 1 to 3. The 
results of experimental measurements of the continuity- equation oscillation (refs. 4 and 
and 5) can be described by equations (7) and (8) of reference 1. This specialized form of 
the continuity equations may be written 


dx _ c „ (v - y) 

dT 71 

(1) 

su-rba-*) 

dr 

(2) 


where, in terms of the parameters used in reference 1, x = N/N q is the dimensionless 



neutral number density, y = N e /N eQ is the dimensionless ratio of electron number den- 

N /N is the ratio of the 
e eo 


sity to its initial (extremal) value, the amplitude index r\ 

average to the extremal electron number density, e = ,/N 0 is the square root of the 

average fraction ionized, and the dimensionless time t = cu o t = 

<ov) ne is the ionization rate parameter for electron- neutral ionization 


V 


N e N 0 ( 0V >ne t ’ where 
The initial con 


ditions of the problem at t = 0 are x(0) = 1.0 and y(0) = 1.0. 

The first terms in the parentheses on the right side of equations (1) and (2) 
follow from the divergence term of the continuity equations when Fick’s law is inappro- 
priate; (ref. 4, appendix B) the second (interaction) terms in the parentheses account 
for the production and loss of each species by the ionization process. 

Approximate solutions to equations (1) and (2) were obtained by Lotka (ref. 6) and 
Volterra (ref. 7) in the linearized limit for which the peak-to-peak amplitude of the fluc- 
tuations of both x and y are small compared with unity. In references 1 to 3 approxi- 
mate closed-form solutions to equations (1) and (2) were obtained in terms of Jacobian 
elliptic functions for the more general case in which only the fluctuations in x(r) were 
assumed small compared to unity. While the method of approximation in references 
1 to 3 may be applied to equations of continuity more general than equations (1) and (2), 
the closed-form solutions it provides to these equations are defective in several re- 
spects. The approximation discussed in references 1 to 3 predicts that oscillatory solu- 
tions to equations (1) and (2) occur for the range 2/3 ^ rj < 1.0, while exact solutions to 
these same equations show that oscillatory solutions exist for all r\ > 0, except the 
steady-state solution at tj = 1. 0. The approximations of reference 1 reproduced the 
exact solutions reasonably well over the range 0. 80 < rj < °°, but a better approximation 
was clearly desirable. 


A NEW APPROXIMATION TO y(x) 

One can obtain an exact relation between x(r) and y(r) by taking the ratio of equa- 
tion (1) to (2), integrating the resulting equation, and applying the initial condition that 
x = 1 when y = 1. This yields a very simple "phase" equation, 

e 2 

(x - 1) - In x = — (1 - y + rj In y) (3) 

V 

Note that a closed curve, indicating oscillating solutions, can be obtained for all r? > 0, 
and that x = 1 at both the maximum and minimum values of y. One of these extremal 
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values of y is the initial value y = 1. If the other extremal value is. denoted y , it is 
apparent from equation (3) that y Q is related to the amplitude index 77 by 


In y„ 

J o 

The parameter y Q is the ratio of the extremal values of y(r). If equation (4) is used to 
replace 77 in equation (1), it becomes evident that if y(y Q ,? 7 )and x(y Q , 77 ) are solutions 
to equations (1) and (2), then so are y(l/y o ,r}/y Q )/y Q and x(l/y Q , ?//y ). As a result of 
this transformation, the solutions for n = 77 > 1 (and y Q > 1) can be obtained from the 
solutions for 77 < 1 (and y Q < 1), and vice versa. The discussion below will be con- 
fined to the case for which 77 < 1 , and y Q = Y min /y max < 1 - 




Dimensionless time, t 

(al rt * 0 . 20 . (b) 17 = 0. 70. 

Figure 1. - Comparison of exact solutions y exac t(T), obtained from the nonlinear pair of continuity equations (eqs. (1) and (2)1, 
with the approximations yj(T) (from eq. (10)), y^h (from eq. (12) for 17 ), and y 4 (t) (from eq. (17) or rj = rj/y 0 ). e ■ 0.01. 


The exact solutions to equations (1) and (2) for y(r) with e = 0. 01 are shown in fig- 
ure 1 for two values of 77 . These solutions were obtained by numerical integration on 
a computer. In figure 2 are plotted the frequency and the average value of the solution 
77 as functions of the peak-to-peak ratio y. These exact solutions will be compared 
with several approximate solutions for y(r) given subsequently. 

If it is assumed that the neutral density fluctuates only very little compared to its 
average value, then x is near 1, and the left side of equation (3) approaches 
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Peak-to-peak amplitude ratio, y m i n /y max 

(a) Average amplitude rj = y/y~ ax . (Amplitudes 0 f (b) Dimensionless period of oscillation. 

y ex act- yp and V 4 are equal. ) 

Figure 2 . - Average amplitude and period of oscillation as functions of peak-to^peak amplitude ratio for the exact and approximate 
solutions. 

o 

1/2 (x - 1) . This requirement of small amplitude oscillations of x(r) can be guaranteed 
by proper choice of e and 77 . The right side of equation (3) takes its maximum value 
at y = rj. Therefore, if one chooses e and rj such that 


2(1 -77-4-77 In rj) 
V 


1/2 


« 1 


( 5 ) 


then one is assured that 


(x - 1 ) - In x « 1 


( 6 ) 
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which implies that x - 1 « 1. After expanding the left side of equation (3), one may 
solve for (x - 1) and obtain 


(x - 1) « ±e 


2 fa In y - y + 1) 


V 


f-i 


= ±e i- ffa, y) 

V 


1/2 


(7) 


Substituting this approximation into equation (2) yields an approximate equation in terms 
of the single variable y(r), 


/ \2 

— (— J « - [v In y - (y - 1 )] = - f fa, y) ( 8 ) 

y2 \dr/ Tj ■ t? 


This approximate equation for y may be integrated by adopting an integrable ap- 
proximation to f(y , y): 


i 1 (y 0 ,y) - -A' V + V 1 ") (in y o - in y)ln y ( 9 ) 

( ln o 2 

which satisfies the requirement that f^ have the same maximum as f, and be zero at 
y = 1 and y = y . The solution obtained by substituting equation (9) into equation ( 8 ) is 

y 1 ( T ) = y< 1 - C 0 S “i T)/2 (io) 


where the parameter is 




8(1 - 77 + r? In 7 ?) 


V 




1/2 


= Wj(r?) 


( 11 ) 


Note that y^ (r) reproduces the peak-to-peak amplitude of the exact solutions and is the 

same function regardless of whether 7? or t? is used. This approximation is plotted in 

figures 1 and 2. It does not appear separately on figure 2(a) because equation (9) was 

chosen so that the relation between t? and y . /y is exact. The approximate fre- 

min 

quency given by equation ( 11 ) differs from the exact value by no more than a few percent 
over the range shown. 
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OTHER APPROXIMATIONS TO y(x) 


The preceding approximation may be compared to the two "small amplitude" ap- 
proximations discussed in reference 1. The approximation using 77 may be written in 
the dimensionless form 


y 2 (r) - 1 + Z 3 SN 2 o) 2 t 

where SN is the Jacobian elliptic sine, and 



/ z X 1 / 2 

u) - ) 


2 \ 677 / 

CO I 

HI 

CSI 

tS 3 

(2 - V) + (2/2)) (77 + 2) 

N 

00 

3) 

111 

1 

1 00 

(2-77)-|/[t7-(2/3)](77 + 2) 


with the elliptic modulus given by 



(12) 


(13) 

(14) 


(15) 


(16) 


This approximation is shown in figure 2, and is seen to fail for r\ < 2/3. A second ap- 
proximation taken from reference 1 is obtained by using 77 = 77 /y in the expressions 


where 


Z„Z„SN w„t 

y 3 (r) = 1 + = — 1 — - 

Z 2 “ Z 3 + Z 3 SN W 3 T 



(17) 


( 18 ) 
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^2 ~ 22(17 ) 


(19) 


Z 3 55 ^3(17) 


( 20 ) 



( 21 ) 


This approximation exists over the entire range of 17 . However, the relation between 
77 and y Q differs from the exact value by a maximum of about 20 percent and the fre- 
quency differs by a maximum of about 15 percent. 

The approximations yg^) and y 3 (r) can be improved by defining an equivalent small 

amplitude index 77 , which is that value of 77 which gives the correct peak-to-peak 
s 

oscillation amplitude. This value of 77 = 77 is then substituted in equation (12) for an 
improved approximation y 4 (T), which is shown in figure 11. 


CONCLUDING REMARKS 

In summary, it has been shown that the approximations given in reference 1 can be 
modified to improve their usefulness and accuracy. A new approximation y^(r) has also 
been developed which uses trigonometric functions and gives the frequency as a function 
of the peak-to-peak amplitude ratio y . This approximation is the most accurate of 
those considered, does not require the use of Jacobian elliptic functions, and provides 
sufficient accuracy for most applications. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 7, 1969, 

129-02. 


REFERENCES 

1. Roth, J. Reece: New Mechanism for Low-Frequency Oscillations in Partially Ionized 

Gases. Phys. Fluids, vol. 10, no. 12, Dec. 1967, pp. 2712-2714. 

2. Roth, J. Reece: Periodic, Small- Amplitude Solutions to the Spatially Uniform Plasma 

Continuity Equations. NASA TN D- 4472, 1968. 


7 



3. Roth, J. Reece: Periodic Small-Amplitude Solutions to Volterra's Problem of Two 

Conflicting Populations and Their Application to the Plasma Continuity Equations. 

J. Math. Phys. , vol. 10, no. 8, Aug. 1969, pp. 1412-1414. 

4. Roth, J. Reece: Experimental Observation of Low-Frequency Oscillations Described 

by the Plasma Continuity Equations. Phys. Fluids, vol. 12, no. 1, Jan. 1969, 

pp. 260-262. 

5. Roth, J. Reece: Experimental Observation of Continuity-Equation Oscillations in 

Slightly Ionized Deuterium, Neon and Helium Gas. Plasma Phys. , vol. 11, no. 9, 
Sept. 1969, pp. 763-777. 

6. Lotka, Alfred J. : Undamped Oscillations Derived from the Law of Mass Action. 

J. Am. Chem. Soc. , vol. 42, no. 8, Aug. 1920, pp. 1595-1599. 

7. Volterra, V. : Variations and Fluctuations of the Number of Individuals in Animal 

Species Living Together. International Council for the Exploration of the Sea. 
Journal du Conseil, vol. 3, 1928, pp. 1-51. 


8 


NASA-Langley, 1970 25 E-4960 



National Aeronautics and Space Administration 
Washington, D. C. 20546 

official business FIRST CLASS MAIL 



POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


POSTMASTER: 


If Undeliverable (Section 15S 
Postal Manual ) Do Not Return 


"The aeronautical and space activities of the United Stales shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technolog)' 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained trom: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


